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ABSTRACT 


In  seismic  discrimination  problems  one  is  interested  in  filtering  out  information 
carried  in  certain  waves  from  noise  carried  in  other  waves  With  the  use  of  arrays  of 
sensors,  this  filtering  can  be  performed  in  space  and  time.  It  is  quite  useful  to 
visualize  the  filtering  in  the  Fourier -conjugate  space  of  wave  number  and  frequency, 
and  this  point  of  view  is  developed  here.  Most  of  the  report  is  devoted  to  an  expression 
of  the  pertinent  facts  about  the  propagation  of  elastic  waves  in  the  earth  in  terms  of  the 
frequency-wave  number  diagram  The  remainder  concerns  signal  and  noise  models 
and  the  interpretation  of  filtering  techniques  in  frequency -wave  number  space 
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I.  INTRODUCTION 


The  successful  development  of  statistical  decision  procedures  for  seismic  data 
processing  depends  directly  on  the  accuracy  with  which  the  signal  and  noise  models  em¬ 
ployed  reflect  the  properties  of  the  actual  seismic  waves.  In  detail  these  models  must 
be  empirically  determined,  but  there  are  certain  general  principles  to  which  they  must 
adhere,  and  a  study  of  these  principles  themselves  provides  valuable  insight  into  the 
problem  of  model  construction  The  general  principles  follow  in  turn  from  the  fact  that 
seismic  signals  and  noise  are  elastic  waves  propagating  in  a  nearly  spherically  sym¬ 
metric,  radially  stratified  medium,  and  that  they  are  produced  and  observed  at  (or  very 
near)  the  surface. 

The  salient  features  of  the  elastic  properties  of  the  earth  are  its  high  degree  of 
spherical  symmetry  and  the  fact  that  the  continuous  variation  of  propagation  speeds 
with  radial  distance  is  interrupted  by  a  number  of  discontinuities,  giving  the  medium  a 
spherically  layered  structure.  The  discontinuity  which  separates  the  mantle  from  the 
core  (at  a  radius  of  0,  548R,  where  R  is  the  mean  earth  radius)  strongly  affects  the 
propagation  of  body  waves,  and  must,  of  course,  be  analyzed  with  due  regard  for  the 
spherical  geometry  involved  However  the  Mohorovi^ic  discontinuity  at  the  base  of  the 
crust  and  various  discontinuities  within  the  crust  lie  very  near  the  surface,  and  they 
give  rise  to  waves  confined  to  the  immediate  vicinity  of  the  surface  Asa  result,  the 
propagation  features  associated  with  these  layers  are  well  described  by  an  equivalent 
plane -layered  model  earth,  consisting  of  one  or  more  layers  overlying  an  infinite  half- 
space. 

A  source  in  the  layered  region  will  send  energy  along  the  surface  in  various  wave¬ 
guide  modes,  some  of  which  (trapped  modes  such  as  Rayleigh  and  Love  waves)  confine 
the  energy  flux  rigorously  to  the  surface,  while  others  (leaking  modes  such  as  PL)  leak 
energy  slowly  down  into  the  half-space  In  addition  to  these  waves,  the  source  will  send 
ordinary  body  waves  down  into  the  half-space.  We  then  interpret  the  half  space  as  the 
interior  of  the  earth,  and  superpose  on  the  plane -layered  model  the  refractive  and  other 
effects  that  will  allow  energy  to  re-emerge  at  the  surface. 

The  theoretical  determination  of  the  elastic  field  produced  by  a  given  source  dis¬ 
tribution  constitutes,  in  general,  a  very  difficult  wave -propagation,  or  diffraction  problem. 


One  must  solve  the  elastic  wave  equation  for  prescribed  sources  and  given  boundary 
conditions  (at  the  earth's  surface  and  at  the  layer  interfaces).  However,  by  means  of  a 
suitable  Green  s  function,  this  solution  for  given  sources  can  be  expressed  as  a  super¬ 
position  of  simpler  solutions,  the  eigenfunctions  of  the  problem  These  eigenfunctions 
are  generalized  plane  waves  (and  they  reduce  to  simple  plane  waves  for  an  infinite 
homogeneous  medium),  with  harmonic  time  dependence,  and  they  satisfy  the  boundary 
conditions  as  well  as  the  source -free  elastic  wave  equation  The  eigenfunctions  con¬ 
stitute  a  three -parameter  family  of  functions.  One  parameter  is  frequency  (which 
governs  the  time -dependence)  and  the  other  two  are  often  taken  to  be  the  plane -wave 
direction  angles.  However,  in  the  seismic  problem  with  its  horizontal  plane  bounding 
surface,  it  is  more  convenient  to  use  azimuth  and  horizontal  wave  number  (defined 
below)  for  parametrization  The  representation  of  an  arbitrary  elastic  field  as  a  super¬ 
position  of  eigenfunctions  is  analogous  to  the  use  of  the  Fourier  integral  for  functions  of 
time.  However,  here  the  analog  of  the  transform  function  of  frequency  is  a  function  of 
the  three  parameters:  frequency,  horizontal  wave  number  and  azimuth  For  remote 
sources  in  a  plane -layered  medium,  the  role  of  the  azimuth  parameter  is  trivial,  and 
it  relates  only  to  the  source  distribution  However  the  frequency,  uu ,  and  horizontal 
wave  number,  K,  are  intimately  related  to  the  properties  of  the  medium  as  well  as  the 
source,  hence  the  representation  of  elastic  waves  in  terms  of  their  density  in  the  uu  —  K 
plane  provides  a  convenient  focal  point  for  the  discussion  of  the  general  principles  men¬ 
tioned  above.  In  the  case  of  seismic  noise,  the  appropriate  description  is  in  terms  of  a 
spectral  density  in  the  uu  —  H  piane. 

In  the  following  sections  we  discuss  the  form  of  the  basic  eigenfunctions  and  the 
significance  of  regions  of  the  uu— k  plane  beginning  with  a  simple  example  and  ending  with 
a  qualitative  discussion  of  the  actual  earth  A  final  section  then  deals  with  models  of 
seismic  signals  and  noise,  emphasizing  the  fundamental  restrictions  imposed  by  the 
physical  nature  of  the  medium  The  discussion  is  expository  throughout,  and  contains 
no  new  results.  Sections  11,  III  and  IV  are  intended  as  aids  in  the  understanding  of 
some  of  the  basic  facts  of  theoretical  seismology.  *  Section  V  provides  the  framework 
within  which  models  can  be  made  and  processing  techniques  can  be  evolved  and  analyzed. 

*  See  References  1  and  2  for  further  details 
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II.  ACOUSTIC  WAVES  IN  AN  INFINITE  HALF -SPACE 


Let  the  scalar  quantity,  cp(r,t),  be  the  velocity  potential  of  an  acoustic  field, 
satisfying  the  wave  equation 


V2cp(r,t)-  \  8  =  0  (2.1) 

cz  ot 

— »  — *  — »  3  cpf  r  t) 

The  velocity  field  itself  is  v(r,t)  =  and  the  pressure  is  p(r,t)  =-p - ^ —  , 

where  p  is  the  density  of  the  medium. 

For  an  infinite  homogeneous  space,  the  basic  eigenfunctions  are  simple  plane 

waves: 


cp(r,t) 


.  i(k  r-'JUt) 

A  e 


(2.2) 


where  A  is  a  complex  constant  and  the  real  and  imaginary  parts  of  (2.2)  are  each  solu¬ 
tions  of  (2.  1)  provided 


|k|2  -  (oj/c)2  =  0 


(2.3) 


According  to  (2.3),  the  eigenfunctions  are  a  three -parameter  family,  which  may  he 
indexed  by  the  three  components  of  k  (in  which  case  uu  =  |k|c),  or  by  uu  and  the  two 
direction  cosines  of  k  (in  which  case  |k|  =  uu/c). 

Although  we  always  keep  the  frequency  real,  Eq*  (2  3)  does  not  force  the  com¬ 
ponents  of  k  to  be  real.  For  example,  kg  (the  z-component  of  k)  might  be  pure  imaginary, 
so  that  kg  is  negative,  without  violating  (2.3)  The  solution  (2.2)  would  then  have  the 
form 


cp(r,t)  =  Aexp(±|k3|z)  exp  [  i(kjX  +  k  y  -  uut)  ] , 

and  would  increase  exponentially  in  either  the  positive  or  negative  z -direction.  Since  this 
cannot  be  permitted  in  an  infinite  medium,  it  can  be  seen  that  only  k-vectors  with  real 
components  are  possible  solutions  of  (2.3).  Such  solutions  are  called  real  plane  waves 
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In  the  problem  of  the  half-space,  we  suppose  that  the  boundary  is  the  plane, z  =  0, 
and  that  the  medium  lies  in  the  upper  space,  z  a  0.  We  now  form  the  eigenfunctions  of 
this  problem  out  of  the  plane  waves  (2.2).  First  we  solve  Eq.  (2.3)  for  k^: 

k  =  ±  7(-)2-h2  (2.4) 

o  c 


where  h,  defined  by 


K 


(2.5) 


is  the  "horizontal  wavenumber  "  We  let  k  be  the  vector,  in  the  x-y  plane,  whose  com¬ 
ponents  are  and  k^,  and  whose  magnitude  is  H.  We  also  let  p  be  the  projection  of  the 
position  vector,  r,  on  the  x-y  plane,  so  that 


kjX  +  k2y 


K  p 


(2.6) 


Then  the  plane  wave  (2.  2)  can  be  written  in  the  form 


cp(r,t) 


A  e 


ik3z  i (h  p  ~ out) 


(2.7) 


which  is  especially  suited  to  a  treatment  winch  emphasizes  the  values  of  the  field  on  the 
boundary  plane  itself.  In  fact,  we  now  choose  to  index  the  eigenfunctions  by  the  param¬ 
eters  frequency  (uu),  horizontal  wave  number  (h)  and  azimuth  (direction  of  h). 

According  to  (2.4)  there  are  two  plane  waves  for  each  set  of  parameter  values. 
However,  only  a  single  combination  of  these  two  will  fit  the  boundary  condition  on  the 
plane  z  =  0.  It  can  be  shown  that  the  field  is  uniquely  determined  by  either  of  two  con¬ 
ditions  at  the  boundary 


or 


a)  cp(r,t)  =  0 


b) 


d^(r,t)  _ 

dz 


at  z  =  0 

at  z  =  0 


(2.8a) 

(2.8b) 
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Sen 

It  follows  from  condition  a)  that  vanishes  at  z  =  0,  hence  the  pressure  is  zero  on 
the  boundary.  This  is  called  a  "soft  boundary"  condition  and  corresponds  to  the  pres¬ 
ence  of  a  free  surface.  The  other  condition,  (2.8b),  says  that  v  ,  the  z -component  of 
velocity,  vanishes  on  the  boundary,  and  this  is  called  a  "hard  boundary"  condition, 
since  the  boundary  surface  is  not  allowed  to  move.  It  is  also  possible  to  determine  a 
unique  solution  by  means  of  a  mixed  boundary  condition 


a<p£.t)  +  b  =  0 


on  z  =  0, 


or  equivalently,  in  terms  of  pressure  and  velocity, 

p(x,y,  0,t)  = 

v^(x, y,  0, t) 

This  is  called  an  "impedance  boundary  condition,  "  and  Z  is  the  mechanical  independence 
(per  unit  area)  of  the  surface.  It  is  useful,  for  example,  when  the  bounding  surface  is 
a  piston  of  finite  mass  and  compliance. 

To  satisfy  any  of  these  boundary  conditions  for  a  fixed  parameter  set  (uu,h),  a 
combination  of  the  two  real  plane  waves  (2.7)  must  be  used.  If  we  take,  for  k^,  the 
positive  root 


k3  =  y  (uu/c)2  -  h2  , 


then  we  obtain  the  basic  eigenfunctions: 


a)  Soft  boundary 


.  .  ikr,z  i(K'p-uut)  .  -iloiz  i(H*  p-uut) 

cp(r,t)  =  A  e  J  e  -  A  e  Je 


=  (2iA)  sin  (k^z)  e 


i(K-  p-uut) 
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b)  Hard  boundary 


cp(r,t)  =  A  e 


=  2A  cos  (k^z)  e 


In  these  equations,  A  is  complex,  and  the  real  and  imaginary  parts  of  cp(r,t)  provide 
independent  eigenfunctions  which  differ  only  in  the  time  origin.  For  example,  in  the 
soft  boundary  case  (which  corresponds  to  the  seismic  problem)  the  two  real  solutions 
are 


sin  (k^z)  cos  (ft-  p-uut) 


(2.9) 


and 


(2.  9') 


The  representation  in  terms  of  the  difference  of  plane  waves  shows  that  each  eigen¬ 
function  may  be  thought  of  as  an  incident  and  reflected  wave,  with  equal  angles  of 
incidence  and  reflection,  and  a  reflection  coefficient  (for  velocity)  of  (—  1).  For 
economy  of  writing,  we  use  only  the  soft  boundary  condition  from  now  on,  although 
perfectly  analogous  results  hold  for  the  hard  boundary  or  (more  generally)  the  impedance 
boundary  condition  problems. 

It  is  interesting  to  derive  the  particle  motion  that  obtains  for  one  of  our  eigen¬ 
functions,  say  (2.  9).  The  component  of  velocity  in  the  x-y  plane  is  given  by 


-  H  sin  (k,^z)  sin  (ft  p-uut), 


which  vanishes  on  the  boundary  The  * 'vertical1 1  component,  however,  is 


which  is  a  maximum  at  the  boundary  Thus  the  particle  motion  at  the  boundary  is 
linearly  polarized,  perpendicular  to  the  boundary,  and  has  the  form  of  a  wave  traveling 
over  the  surface  with  frequency  uu,  apparent  wave  number  n  (wave  length  2tt/k.), 
direction  ft,  and  horizontal  phase  velocity  uu/K,  Since  we  always  have 

k  ^  uu/c  (2.  10) 

for  real  plane  waves  (real  k-vector),  the  horizontal  phase  velocity  is  never  less  than 
the  speed  of  sound 

For  each  point  in  the  uu-H  space  satisfying  (2  10),  we  have  an  eigenfunction 
formed  from  real  plane  waves,.  This  region  clearly  forms  a  cone  about  the  uu-axis,  and 
in  the  uu-H  plane  has  the  form  of  the  shaded  region  in  Fig  1.  In  this  figure,  the  slope  of 
the  boundary  line  is  just  c,  the  sound  speed  in  the  medium  The  shaded  region,  which 
contains  all  eigenfunctions,  is  called  the  continuum . 

The  forbidden  region  in  Fig,  1  is  defined  by  the  condition  K  >  uu/c.  By  Eq.  (2.4), 
this  corresponds  to  plane  waves  with  pure  imaginary  values  of  We  have  seen  that 
these  so-called  complex  plane  waves  were  excluded  in  infinite  space  because  they 
become  infinitely  large  in  the  ±  z  direction  However,  for  the  half -space,  we  always 
have  z  ^  0,  so  that  the  case  k  =  ip  where 

o 

p  =  +/k2  -  (oj/c)2  (2.11) 

is  not  excluded.  For  this  case, 

cp(r,t)  =  e  cos  (k- o-uut)  (2.12) 

and  it  appears  that  we  have  another  solution.  However,  without  the  other  complex  plane 

-f-pz 

wave  (which  has  z -dependence  e  ),  we  cannot  fit  any  of  the  boundary  conditions  dis¬ 
cussed  above.  Nevertheless,  because  of  the  importance  of  such  waves  in  the  elastic 
case,  it  is  interesting  to  derive  the  associated  particle  motion.  We  take  the  x-axis 
along  h,  so  that  the  disturbance  propagates  along  the  x-axis  with  phase  velocity  (uu/h)  <  c. 
Then 


7 


sin  (Hx-oit) 


dy(r,t) 
x  9x 

and 

Scp(r,t) 

z  dz 

These  components  of  velocity  are  out  of  phase,  so  that  the  resultant,  although  a  pure 
sound  wave,  is  not  longitudinal  in  the  ordinary  sense.  The  fact  that  the  z -component 
of  velocity  is  out  of  phase  with  the  pressure, 

/  v  dcp(r,t)  -pz  .  .  . 

p(r,t)  =  —  p  — ^ -  =  puue  sin  (Hx-tut)  , 

implies  no  flux  of  energy  in  the  z -direction.  At  a  given  point  on  the  boundary,  the 
particle  trajectory  is  given  by  the  parametric  equations 

x(t)  =  —  H  cos  U)(t-to) 

(2.  13) 

z(t)  =  —  psinuu(t-to) 

The  particle  motion  is  elliptical  (prograde),  with  major  and  minor  ellipse  axes  along  the 
x-  and  z-axes. 


K  e 


-pz 


p  e  cos  (Hx-tut) 
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III.  ELASTIC  WAVES  IN  AN  INFINITE  HALF -SPACE 


Elastic  waves  are  described  by  a  vector  field,  u(r,t), which  represents  the  dis¬ 
placement  from  equilibrium  due  to  the  elastic  disturbance.  Instead  of  discussing  the 
rather  complicated  equation  which  u  obeys,  we  make  use  of  that  fact  that  any  elastic 
field  may  be  expressed  as  the  sum  of  two  fields: 


u(r ,  t)  =  u  (r ,  t)  +  u„(r ,  t) 
P  & 


each  of  which  satisfies  a  simple  wave  equation: 


2  -  .  Id  up(r,  t)  _ 

v  u  (r,t)  — o  — Irr  —  =  0 
P  az  dtz 


V  u  (r,t)  - 


J_  d  us(r,t)  _ 


|3Z  dt^ 


=  0 


and  an  auxiliary  condition: 


V  x  U  (r,t)  =  0 
P 


V  u  (r,t)  =  0 


(3.1) 


(3.2a) 


(3.2b) 


(3.3) 

(3-4) 


The  partial  field  u  (r,t)  is  irrotational ,  describes  dilatational  disturbances,  and  prop- 
P 

agates  with  phase  velocity  a  Plane  waves  of  this  type  are  called  P-waves.  The  partial 
field  u  (r,t)  is  solenoidal,  describes  equi-voluminal  disturbances  (pure  shear),  and 
propagates  with  phase  velocity  0  <  a  Plane  waves  of  this  type  are  called  S -waves 

It  follows  from  (3.  3)  that  u  (r,t)  can  be  derived  from  a  scalar  potential  function, 

'J'(r.t) 


Up(r ,  t)  =  V'lKr.t)  , 


which  satisfies  the  wave  equation  with  velocity  a-  The  time -derivative , 
entirely  equivalent  to  the  velocity -potential ,  cp(r,t),  of  Section  II.  The  u 
identical  in  properties  with  the  acoustic  velocity  field  discussed  above, 
the  basic  Up-  waves  are  real  plane  waves: 


d'l>(r,t) 

dt 

field  is  therefore 


P 

In  infinite  space, 
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(3.5) 


u  (r,t)  =  U  e 
P  P 


i(lc  r-uut) 


where 


|k|2  -  (uu/a)2  =  0 


(3.6) 


and  the  components  of  k  are  real.  For  given  values  of  uu  and  x,  there  are  two  waves, 

— *  i 

(u  ),  one  corresponding  to  each  sign  of  kg.  According  to  condition  (3.3),  which  must 
hold  everywhere  and  for  all  time,  the  polarization  vector ,  U^,  satisfies 


k  X  U  =0 
P 


(3.7) 


Although  may  be  complex,  it  follows  that  the  real  and  imaginary  parts  are  each  par 

allel  to  k,  so  that  U  itself  has  the  form 
P 


U  =  A  k 
P 


(3.8) 


where  A  is  complex  (by  redefining  the  time  origin,  we  can  make  A  real).  Equation  (3.8) 
expresses  the  fact  that  P-waves  are  longitudinal. 

In  the  half-space  in  addition  to  real  plane  waves,  with  x  <  uu/a,  we  again  consider 
the  complex  plane  wave,  (u^)°,  which  is  obtained  if  x  >  uu/a  We  can  still  write  this 
wave  in  the  form  (3.  5)  if  we  interpret  k  as  a  complex  vector- 


k  =  x  +  ip  e. 


(3.9) 


where  ez  is  a  unit  vector  in  the  positive  z-direction,  x  is  the  horizontal  component  of  k, 
as  before,  and 


P  =  +  a/h2  -  (uu/a)2 


(3.  10) 
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-  pz 

The  z -dependence  of  this  wave  is  given  by  the  factor  e  *  Now  condition  (3.3)  is  still 
equivalent  to  (3.7)  (although  k  and  are  both  complex),  and  it  may  be  shown  that  (3.7) 
still  implies  (3.8).  In  the  algebraic  sense,  the  wave  is  still  longitudinal,  although  the 
particle  motion  is  quite  different  from  that  of  a  real  P-wave.  Substituting,  we  find  that 


*kei(K'“t) 


A(k  +  ip  ez)  e 


-pt  i(K-p-ujt) 


(3.11) 


We  choose  the  time  origin  to  make  A  real,  set  it  equal  to  unity,  and  take  the  real  part 
of  u.  The  result  is  the  wave 


H  e 


pz 


COS  (H  '  p  —  U)t)  — 


P  ez  e 


-pz 


sin  (k  •  p  —  tut) 


(3.  11') 


which  describes  prograde  elliptical  particle  motion.  In  fact,  except  for  a  constant  fac¬ 
tor,  the  particle  motion  described  by  (3.  IT)  is  identical  with  that  obtained  in  Section  II 
for  complex  acoustic  waves. 

In  infinite  space,  the  basic  ug-  waves  are  real  plane  waves: 


=  Uc 


i(k*  r -out) 


(3.  12) 


where 

|k|2  -  (uu/3)2  =  0 


(3. 13) 


and  the  components  of  k  are  real.  There  are  two  waves,  (u  )^  for  each  set  (uu,h)  and 
for  each  possible  polarization.  Condition  (3.4)  now  reduces  to 


us  =  0 


(3.  14) 


which  means  that  Us  is  an  arbitrary  complex  combination  of  two  independent  real  vectors, 
each  perpendicular  to  k.  Thus  S -waves  have  transverse  particle  motion.  In  the  half- 
space,  it  is  natural  to  choose  these  two  directions  so  that  one  is  perpendicular  to  both 
k  and  ez,  while  the  other  is  in  the  (k  —  ez)  plane  For  simplicity,  we  take  the  x-axis 
along  K,  and  make  use  of  the  unit  vectors  ex  and  e  ,  along  with  x  and  y  axes,  respectively. 
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Then  the  first  type  of  S-wave  is  called  an  SH-wave,  is  polarized  in  the  horizontal  plane, 
and  has  the  form 


sh 


=  A  e 


i(k-r-uut) 


(3.  15) 


for  either  sign  of  k,^  The  other  wave  is  called  an  SV-wave,  is  polarized  in  the  vertical 
plane  containing  k,  and  can  be  written  in  the  form 


u 

sv 


A(k  x  e  )  e 


i(k-  r-uut) 


(3.16) 


again  for  either  sign  of  k^. 

Now  in  the  half-space,  in  addition  to  linearly  polarized  real  plane  waves,  we 
have  the  new  possibility  of  complex  S -waves,  (u^)0  and  (ugv)°.  We  must  have  H  >  uu/3 
(which  exceeds  uu/ct),  and  we  still  use  (3.  12)  with  the  complex  k-vector: 


where 


k  =  H  +  iq  e^ 


n,  c  T  n 
X 


q  =  +  a/x.2  -  (uu/0)2 


(3.17) 


(3.  18) 


It  is  not  difficult  to  show  that  condition  (3.4)  is  still  equivalent  to  (3  14)  which  is  still 
satisfied  for  two  types  of  wave,  again  of  the  form  (3.  15)  and  (3.  16)  (although  there  is  only 
one  complex  SH-wave  and  one  complex  SV-wave,  since  the  sign  of  q  is  fixed).  The  com¬ 
plex  SH-wave  is  still  truly  transverse,  although  the  complex  SV-wave  is  only  formally 
transverse,  in  the  sense  that  (3.  14)  remains  true.  The  complex  SH-wave  has  the 
specific  real  form 

(ugh)°  =  ey  e  qZ  cos  (k  p  —  uut)  (3.19) 

with  simple  harmonic  linearly  polarized  particle  motion  along  the  y-axis.  The  corre¬ 
sponding  form  for  the  complex  SV-wave  is  obtained  by  substituting  (3.  17)  into  (3.  16)  and 
taking  the  real  part  (A  may  be  taken  real  and  ignored): 
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(3.20) 


(usv)°  =  K  ez  e  cos  (Kx  "^)t)  +  qexe^Z  sin  (Kx  -  uut) 

It  may  be  seen  that  (3.  20),  like  (3  IT),  describes  prograde  elliptical  particle  motion 
(had  we  taken  A  to  be  pure  imaginary- -a  shift  of  time  origin  --  Eq.  (30)  would  have 
come  out  resembling  (3.  11')  even  more  closely. )  Thus  complex  P-waves  and  complex 
SV -waves  are  very  similar  and  can  be  combined  without  destroying  the  character  of  the 
particle  motion. 

In  infinite  space,  only  the  real  plane  waves  can  exist  (as  we  have  seen,  there  are 
three  pairs  of  waves).  In  the  half-space  problem,  we  can  expect  the  boundary  condition 
at  z  =  0  to  determine  the  waves,  or  combination  of  waves,  which  can  exist.  The  elastic 
soft -boundary  condition  is  that  the  force  exerted  by  the  medium  on  the  boundary  surface 
(the  "surface  traction1’)  is  zero.  This  force  is  a  vector  which,  for  our  plane  boundary 
at  z  =  0,  is  given  in  terms  of  the  displacement  field  by  the  relation 

f(p,t)  =  X  ez  (v  •  U)  +  [X  +  p  yuz  (3,21) 

The  derivatives  of  u  are  evaluated  at  z  =  0,  uz  is  the  z -component  of  u,  and  the  constants 
\  and  |a  are  the  Lame  coefficients  which  characterize  the  stress-strain  relations  of  the 
medium.  * 

For  SH -waves,  whether  real  or  complex,  y  u  =  0  and  uz  =  0,  so  that  f  is  given 
by  \j.  du/dz  and  has  only  a  y-component.  Thus  the  boundary  condition,  f  =  0,  reduces  to 
a  single  equation  in  this  case  (P-  and  SV -waves  with  the  same  value  of  H  do  not  give  rise 
to  surface  forces  with  y-components).  If  K  <  uu/0,  we  have  two  real  SV -waves:  (usj1)± 
and  the  combination 


+ 


u  =  (ush}  + 


shr 


*  In  terms  of  X,  p,  and  the  density,  p,  the  propagation  speeds  are  given  by 
a2  =  (X.  +  2|a)/p,  p2  =  |_t/p  . 
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results  in  f  =  0.  Thus  a  single  real  SH-wave  can  exist  in  the  half-space,  consisting 
of  an  incident  wave,  (u^)  ,  and  a  reflected  wave,  (u^)  »  with  reflection  coefficient 
equal  to  (+1).  If  H  >  uu/3,  there  is  only  one  complex  SH-wave,  for  which  f  ^  0,  hence 
this  wave  cannot  exist  Thus  the  uu-H  diagram  for  SH-waves  looks  just  like  that  for 
acoustic  waves  (Fig.  1)  except  that  the  dividing  line  has  a  slope  equal  to  the  shear 
velocity,  3. 

We  now  consider  P-  and  SV -waves  together,  since  it  turns  out  that  the  boundary 
conditions  cannot  be  satisfied  with  either  type  alone.  Since  we  are  dealing  with  two 
phase  velocities,  a  and  9,  there  are  three  (uu-h)  regions  to  be  considered,  defined  by 
the  relations 


X  < 

uj./a 

(3.22a) 

uj/a 

on 

3 

V 

V 

(3.22b) 

uj/3 

<  H 

(3.22c) 

These  are  shown  on  Fig  2  The  slope  of  the  line  between  A  and  B  is  a;  the  other  line 
has  slope  6.  In  all  cases,  the  force,  f,  has  two  components  (in  the  x  and  z  directions), 
hence  the  boundary  condition,  f  =  0,  provides  two  homogeneous  equations  to  be  satisfied 
by  the  various  wave  amplitudes 

In  region  A,  both  P-  and  SV -waves  are  real,  hence  there  are  four  waves  in  all. 
The  general  combination 


can  be  made  to  satisfy  the  two  boundary  equations  and  at  the  same  time,  two  of  the  four 
coefficients  can  be  chosen  arbitrarily  For  example,  the  choice  A  =  1,  B  =  0  corresponds 
to  an  incident  P-wave  and  reflected  P-  and  SV-waves.  The  independent  combination 
C  =  1,  D  =  0  corresponds  to  an  incident  SV-wave,  together  with  both  kinds  of  reflected 
waves.  The  reflection  coefficients  for  the  waves  of  the  incident  type  no  longer  have 
magnitude  unity,  because  energy  is  removed  by  mode  conversion  to  the  other  type. 
However,  in  each  of  the  two  cases  given  above,  the  amplitudes  of  both  reflected  waves 
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are  real.  This  results  in  a  linearly  polarized  particle  motion  in  the  vertical  plane. 
Region  A  corresponds  to  all  possible  angles  of  incidence  for  P-waves,  and  the  reflected 
SV-wave  travels  in  a  direction  closer  to  the  normal  to  the  surface,  as  in  Fig.  3.  How¬ 
ever,  for  incident  SV -waves,  region  A  corresponds  to  angles  of  incidence  no  greater 
than  0  ,  the  critical  angle  at  which  the  reflected  P-wave  is  parallel  to  the  boundary;  see 
Fig.  4. 

In  region  B  of  the  uu-h  plane,  there  are  two  real  SV -waves,  but  only  a  single, 
complex,  P-wave.  The  general  case  is  a  superposition' 

u  =  A  (Up)°  +  B(ugv)+  +  C(ugv) 

Again,  the  boundary  equations  can  be  satisfied  and  at  the  same  time  one  of  the  three 
amplitudes  can  be  chosen  arbitrarily.  For  example,  we  may  take  B  =  1,  and  we  have  an 
incident  SV-wave  (with  angle  of  incidence  greater  than  0c),  a  real  reflected  SV-wave 
(angle  of  reflection  equal  to  angle  of  incidence),  and  a  "reflected1*  complex  P-wave. 

In  this  case  the  reflection  coefficients,  A  and  C,  are  complex,  although  |c|  =  1.  The 
particle  motion  is  elliptic  in  the  vertical  plane. 

In  region  C,  we  have  only  two  waves  a  complex  P  and  a  complex  SV  ,  The  com  - 
bination 


u  =  A(up)°  +  B(usv)° 

describes  a  complex  wave  with  elliptical  particle  motion  in  the  vertical  plane.  The 
boundary  equations  are  a  pair  of  homogeneous  equations  in  the  two  amplitudes,  A  and  B. 
Thus  a  non-vanishing  solutions  exists  only  if  the  system  determinant  vanishes,  and  this 
amounts  to  a  relation  between  uu  and  h.  In  fact,  this  relation  simply  gives  the  ratio, 
c  ~  uu/h,  as  a  function  of  the  phase  velocities  a  and  3-  In  other  words,  surface  waves 
of  this  type,  called  Rayleigh  waves,  can  exist  only  for  one  value  of  K  at  each  frequency. 
For  the  homogeneous  half -space  under  study  here,  this  * 'dispersion  relation*'  giving  K 
as  a  function  of  uu,  is  linear0 

K  =  UJ/CR, 
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hence,  Rayleigh  waves  have  a  characteristic  velocity  and  there  is  actually  no  dispersion. 
For  all  values  of  a  and  8,  c^  <  8,  and  for  most  practical  cases,  c^  is  nearly  equal  to 

B. 

We  can  now  summarize  the  state  of  affairs  in  terms  of  two  uu-H  diagrams,  one 
for  SH -waves  and  one  for  both  P-  and  SV -waves.  These  are  shown  in  Figs.  5  and  6,  in 
which  the  labels  on  the  lines  give  their  slopes.  For  all  points  in  the  uu-H  plane  to  the 
left  of  the  line  S,  in  Fig.  5,  we  obtain  an  eigenfunction  representing  real  SH-waves.  In 
the  same  region  of  Fig.  6,  we  get  eigenfunctions  representing  combinations  of  P-  and 
SV -waves,  with  at  least  two  waves  real.  To  the  right  of  line  3,  in  Fig.  6,  we  get  a 
family  of  eigenfunctions,  all  of  which  correspond  to  points  on  the  line  c^,  representing 
Rayleigh  waves,  made  up  of  complex  P  and  SV -waves.  All  SH-waves  are  linearly 
polarized  (transversely)  in  the  horizontal  plane.  Eigenfunctions  of  the  P-SV  type  are 
linearly  polarized  in  the  vertical  plane  for  points  to  the  left  of  line  a.  To  the  right  of 
line  a,  the  particle  motion  is  elliptic*  in  the  vertical  plane.  For  any  eigenfunction  in 
the  uu-H  plane,  any  azimuth  is  possible .  This  exhausts  the  possibilities  for  the  case  of 
a  homogeneous  elastic  half- space. 


*  In  Rayleigh  waves,  the  P  and  SV  components  are  separately  prograde  in  partical 
motion  at  the  surface,  although  the  combination  may  be  retrograde  (as  in  the  actual 
earth). 
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IV,  ELASTIC  WAVES  IN  A  LAYERED  HALF-SPACE 


As  mentioned  in  the  Introduction,  most  of  the  observed  features  of  seismic  wave 
propagation  can  be  accounted  for  by  a  hybrid  treatment  which  consists  of  ray  methods 
for  body -wave  phases  traveling  through  the  interior  of  the  earth,  and  the  rigorous  dif¬ 
fraction  theory  of  plane -layered  media  for  surface  waves.  The  methods  of  the  previous 
section  apply  equally  well  to  the  case  of  n  plane  layers  (of  arbitrary  type  and  thickness) 
overlying  a  half-space,  although  the  analysis  is  tedious  and  the  results  are  complicated 
in  detail.  Many  of  the  observed  features  of  surface  waves  can  be  explained  by  assuming 
only  one,  or  at  most  two,  layers  With  infinite  parallel  plane  layers,  moreover,  SH- 
waves  never  mix  with  P-SV -waves,  hence  the  two  cases  can  be  handled  separately,  and 
two  uu-k  diagrams  still  suffice  to  describe  the  possible  eigenfunctions.  Since  the  surface 
wave  or  body-wave  character  of  an  eigenfunction  depends  entirely  on  its  behavior  in  the 
half-space,  one  need  only  plot  the  line  uu  =  0k  on  the  SH  diagram,  and  the  two  lines  uu  =  ax 
and  uu  =  3k  on  the  P-SV  diagram,  to  separate  these  two  classes  of  waves  Here,  a  and  (3 
are  the  P-wave  and  S-wave  velocities,  respectively,  in  the  underlying  infinite  half  space 

First,  we  discuss  SH -waves.  If  k  <  uu/0,  we  are  in  the  continuum,  and  real  plane 
SH -waves  can  exist  in  the  underlying  half  space  Moreover,  there  are  two  of  these 
waves,  one  traveling  toward  and  one  away  from  the  surface  *  In  each  layer  (of  finite 
thickness)  two  waves  can  exist  as  well.  Depending  on  the  values  of  uu  and  K  (and  the 
properties  of  the  layer),  the  waves  in  a  given  layer  will  either  be  real  or  complex* 

Since  ail  SH  displacements  and  forces  are  in  a  single  direction  (for  fixed  uu  and  H),the 
boundary  conditions  at  an  interface  (continuity  of  displacement  and  surface  force)  amount 
to  just  two  equations  Again,  at  the  free  surface,  z  =  0,  the  boundary  condition  yields 
one  equation.  Thus,  for  the  case  of  n  layers  over  a  half  space,  there  are  2n  +  2  wave 
amplitudes  and  a  total  of  2n  T  1  homogeneous  boundary  equations.  Thus  we  can  always 
solve  the  system,  while  specifying  one  amplitude  (say  that  of  the  incident  wave  in  the 
half-space)  arbitrarily.  The  solution  describes  the  multiple  reflection  of  an  incident 
real  SH-wave  by  the  set  of  layers 

Now,  if  k  >  uu/B;  we  have  a  single  permitted  complex  SH-wave  in  the  half-space 
and  two  waves  in  each  layer  The  boundary  conditions  reduce  to  2n  f  1  homogeneous 
equations  in  2n  +  1  unknowns,  and  solutions  are  possible  only  if  the  determinant  of 


17 


this  system  of  equations  vanishes.  This  condition  is  called  the  "period  equation";  its 
solution  specifies  mas  a  multivalued  function  of  h.  Thus  we  have  a  series  of  waveguide 
modes,  called  Love  waves,  each  characterized  by  a  dispersion  relation  (solution  of  the 
"period  equation")  giving  uu  as  a  function  of  H.  From  this  relation  we  can  derive  the 
equations  for  phase  velocity  (uu/h)  and  group  velocity  (duu/dH)  as  functions  of  H. 

For  the  case  of  a  single  layer,  the  behavior  of  Love  waves  is  fairly  simple.  They 
can  exist  only  if  the  shear  velocity  (3)  in  the  half-space  exceeds  that  (S')  in  the  layer. 
Also,  for  a  given  frequency,  uu,  the  system  determinant  has  zeros  only  for  values  of  H  in 
the  range 


UJ 

3 


K 


uj 

3’ 


(4.1) 


Thus,  on  the  SH  uu-H  diagram,  the  dispersion  curves,  which  also  represent  the  loci  of 
surface-wave  eigenfunctions,  lie  entirely  between  the  lines  labelled  3  and  3'  in  Fig.  7. 

In  this  wedge-shaped  region  the  waves  in  the  layer  are  real,  while  those  in  the  half -space 
are  complex.  The  period  equation  admits  an  infinite  sequence  of  modes  (each  corre¬ 
sponding  to  an  integral  number  of  nodal  planes  in  the  layer),  each  confined  to  a  dispersion 
curve  of  uu  versus  H,  as  shown  in  Fig.  7.  The  first  mode  exists  for  all  frequencies,  and 
the  higher  modes  exist  for  frequencies  above  the  cut-off  frequencies: 

fn  *  (8/H)  t(3/3’)2-  1]'1/2  (4.2) 


where  H  is  the  thickness  of  the  layer.  Each  mode  propagates  with  phase  velocity  3  near 
cut-off,  and  with  phase  velocity  3'  as  uu  becomes  infinite.  In  other  words,  the  phase 
velocity  is  that  of  the  layer  for  wavelengths  small  compared  to  H,  and  that  of  the  half¬ 
space  in  the  opposite  case. 

The  behavior  of  P—SV -waves  is  more  complicated  in  detail.  We  now  have  four 
waves  in  each  layer  and  four  continuity  equations  at  each  interface.  Again  there  are 
solutions  in  the  continuum  portion  of  the  uu-H  diagram,  representing  incident  P-  or  SV- 
waves  and  reflected  waves  of  both  types.  There  is  also  the  region  where  P-waves  are 
complex  and  SV-waves  real  in  the  half-space,  as  before.  For  values  of  h  greater  than 
uu/3,  both  waves  in  the  half -space  are  complex,  and  the  boundary  conditions  for  n  layers 
reduce  to  4n  +  2  homogeneous  equations  in  4n  +  2  unknown  wave  amplitudes.  The  resulting 
period  equation  again  yields  an  infinite  series  of  modes,  each  with  its  characteristic 
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equation  of  uu  versus  H  For  a  single  layer,  we  find  modes  entering  in  pairs  as  uu  in- 
ereases  past  successive  cut-off  values,  and  all  the  dispersion  relations  are  curves  in 
the  uu-H  plane.  Hence  Rayleigh -wave  propagation  is  now  dispersive,  like  Love -wave 
propagation,  so  that  the  non -dispersive  Rayleigh  waves  of  the  pure  half-space  appear  as 
a  very  special  case.  Typical  curves  are  sketched  qualitatively  in  Fig.  8  for  P—SV -waves 
in  a  medium  consisting  of  a  single  layer  over  a  half-space.  The  three  straight  lines 
are  labelled  by  their  slopes,  a,  6,  and  cD ,  equal  to  the  phase  velocities  of  P,  S  and 
Rayleigh  waves ,  respectively,  in  the  half-space.  (The  indicated  asymptotic  slopes  are 
S',  the  shear  velocity  in  the  layer,  and  c^’,  the  Rayleigh  wave  velocity  in  the  layer.) 

Up  to  this  point,  we  have  considered  only  eigenfunctions  made  up  of  plane  wave 
functions  proportional  to 

i(K-  p— cut) 
e 


for  real  values  of  h  and  uu.  The  z  dependence  is  then  either  trigonometric  or  exponential, 
depending  on  the  relative  magnitudes  of  h  and  uu.  No  other  solutions  are  possible  which 
are  finite  throughout  all  space  (or  half-space)  and  for  all  time.  However,  the  wave 
equation  (for  P  waves,  for  example)  is  satisfied  by  the  plane  wave 


up(r,t)  =  Up  e 


i(k3Z  +  K‘  p— uut) 


(4.3) 


for  complex  values  of  k^,  H,  and  uu,  so  long  as 

2  2  2 
k3  +  K  -  (uu/a)  -  0 


(4  4) 


For  the  layered  half-space  problem  it  is  found  that  the  boundary  conditions  can  still  be 

- 1 

satisfied  for  combinations  of  waves  with  real  uu  and  complex  h  To  put  it  another  way, 

the  same  conditions  (vanishing  determinant)  which  yield  the  dispersion  relations  for 

— < 

Love  and  Rayleigh  waves  admit  solutions  for  real  uu  and  complex  h 

In  order  to  see  the  physical  nature  of  such  waves,  suppose  that  h  is  in  the  x 
direction  and  that 


h  =  a  +  ib 
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where  a  and  b  are  real.  Then  (unless  a  =  0)  H  is  complex,  and  (4.4)  implies  that  k  , 

o 

and  hence  k  itself,  is  complex.  If  we  put 
o 

=  c  +  id 


then  the  plane  wave  (4.3)  becomes 

-*  r*  x  rt  —  dz  +  icz  —  bx  +  lax  —  iuut  ..  rx 

up(r,t)  =  Up  e  (4’5) 

Formula  (4.  5)  describes  a  complex  plane  wave  of  a  more  general  type  than  those  con¬ 
sidered  before.  The  disturbance  propagates,  with  attenuation,  in  both  x  and  z  directions. 
It  turns  out  that  solutions  exist  for  positive  values  of  a,  b,  and  c,  and  negative  d.  Thus 
energy  travels  along  the  surface  (positive  x  direction)  with  exponentially  decreasing 
amplitude,  while  energy  is  also  propagating  into  the  half-space  (positive  z  direction) 
with  exponentially  increasing  amplitude.  These  waves  are  not  true  surface  waves 
because  of  the  positive  flux  of  energy  out  of  the  surface -layer  waveguide  into  the  half- 
space,  and  for  this  reason  are  called  leaky  modes.  Such  waves  obviously  cannot  exist 
in  the  steady  state  (real  uu)  because  they  have  unbounded  amplitude.  However,  they  have 
a  definite  place  in  transient  wave  propagation  phenomena  in  which  they  enter  the  general 
eigenfunction  expansion  with  complex  values  of  u).  These  modes  may  be  indicated  on  the 
uu-H  diagram  (for  real  uu)  by  plotting  uu  versus  the  real  part  of  H.  They  appear  as  a  series 
of  dispersion  curves  located  in  the  continuum  portion  of  the  diagram. 
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V.  SIGNAL  AND  NOISE  MODELS 


Let  us  suppose,  at  first,  that  we  require  signal  and  noise  models  suitable  for 
the  study  of  a  surface  array  of  vertical -component  sensors.  The  vertical  component 
of  earth  displacement,  uz,  as  a  function  of  position  on  the  surface,  p,  and  time,  t,  can 
be  written  in  the  form 

uz(p,t)  =  JJJ  i|i(w,K)  e1^  P  ^dHduu  (5.1) 

—  00 

Formula  (5.  1)  is  introduced  here  simply  as  a  three-dimensional  Fourier  transform, 
where  frequency  (uu)  is  conjugate  to  time,  and  the  components  of  horizontal  wave  number 
(h)  are  conjugate  to  the  horizontal  position  coordinates.  Since  uz  is  a  real  quantity,  it 
follows  that 


♦  (-«),  -  H)  =  t(uo,K)  (5.2) 

where  the  bar  denotes  complex  conjugate.  As  a  result  of  (5.2)  we  can  rewrite  (5.  1)  in 
the  form 

00  00 

u  (p,t)  =  Jj  dHj  dju  {  X'(iu,K)  cos  (k  •  p— uu t)  -  Y'(uu,  h)  sin  (h  •  p-uut)  }  (5.3) 

—  CO  0 

where  X'  and  Y'  are  real,  and 

X’(u),K)  +  i  Y’(uj,H)  =  2i|r(uu,K)  (5.4) 

If  we  use  polar  coordinates  in  the  H-plane,  so  that  H  is  the  length  of  K  and  a  is  the  polar 
angle,  and  also  polar  coordinates  in  space,  with  p  =  |p|  and  azimuth  angle  0,  then  (5.3) 
can  be  written 

2TT  CO 

uz(p,0,t)  =  J  d0  JJduudH  {X(UU,H,a)  COS  (Hp  COS  (a-0)  -  out] 

0  0 

—  Y(u),H,a)  sin  [Kp  COS  (a-0)  -  tut]  }  (5.5) 

where  X(uu,H,a)  =  kX'(uu,h)  and  Y(uu,K,a)  =  kY'(uj,h). 


21 


Now  that  the  field  is  represented  as  an  integral  over  azimuth  (a)  in  the  first 
quadrant  of  the  uu-H  plane,  we  may  interpret  X  and  Y  as  amplitudes  in  the  same  uu-H 
space  discussed  in  earlier  sections.  In  other  words,  form  (5.5)  makes  it  possible  to 
inject  into  the  model  such  physical  assumptions  as  "signal  is  P-wave  only"  or  "noise 
is  all  in  the  first  Rayleigh  mode  ’  etc.  ,  simply  by  limiting  non-zero  values  of  X  and  Y 
to  appropriate  portions  of  the  uu-H  space.  Generally,  signals  contain  many  phases,  and 
noise  is  present  in  several  surface -wave  modes,  hence  adequate  models  will  require 
amplitudes  (X  and  Y)  which  are  sums  of  several  terms,  each  relatively  sharply  con¬ 
fined  to  a  portion  of  the  uu-H  space. 

The  representation  (5.5)  is  directly  suitable  to  model  deterministic  signals  of 
known  or  unknown  form.  However,  to  describe  noise,  or  to  derive  a  random  model  for 
signals,  we  begin  with  a  statistical  specification  and  proceed  to  derive  an  analog  of 
(5.5).  Let  us  assume  that  uz  is  a  random  function  of  p  and  t  with  mean  zero: 

Euz(p,t)  -  0  (5.6) 

The  symbol  E  stands  for  expectation,  or  ensemble  average.  We  also  assume  that  the 
noise  is  (wide -sense)  stationary  in  both  space  and  time,  so  that  its  covariance  function, 
T,  has  the  form 

E  uz(p,t)  uz(p’,t’)  =  T(p-p',  t-t')  (5.7) 

(The  assumption  of  stationarity  is  a  convenience,  not  a  necessity;  very  similar  results 
can  be  obtained  without  its  use,  but  the  analysis  is  more  complex. )  The  covariance 
function,  T(p,t),  always  has  a  spectrum;  we  assume  it  has  a  spectral  density,  G'(uu,h) 

r<p,t)  =  /Jj  G'ta.K)  (5.8) 

—  00 

As  before,  this  relation  is  a  purely  mathematical  one,  essentially  a  Fourier  transform, 
and  we  choose  to  call  the  variables  conjugate  to  p  and  t  by  the  names  h  and  0),  respec¬ 
tively. 
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Because  T(p,t)  is  a  covariance,  G'(uj,h)  is  real  and  non-negative.  Because  uz, 
and  hence  r(p,t)  is  real,  we  have 


G'(—  uj,  —  h)  —  G'(uu,  k) 


(5  9) 


This  symmetry  allows  us  to  write  (5.8)  as  an  integral  over  positive  frequencies  only: 

00  00 

r(p,t)  =  2j  duu  Jj  dH  G’(uu,  h)  cos  (k  •  p  —  uut)  (5.10) 

0 

Finally,  in  polar  coordinates,  we  put 


G((ju,H,a)  =  2kG’(uu,k) 

and  find 

2tt  °° 

r(p,0,  t)  =  \  da  J  duu  dx  G(uu,  k,  a)  cos  [  Kp  cos  (a-0)  —  uut]  (5.11) 

0  0 


The  well-known  spectral  representation  theorem  for  stationary  random  functions 
of  time  has  a  direct  analog  for  stationary  random  functions  of  several  variables.  When 
applied  to  the  random  function  uz(p,t)  and  expressed  in  polar  coordinates,  the  represen¬ 
tation  is 


2  TT  00 

uz(p,9,t)  =  J  {cos  [hdcos  (a-0)  —  tut]  d5(aj,K,a) 

0  0 

—  sin  [Kp  cos  (a-0)  -  uut]  dr|(cij,x,a)  }  (5.  12) 

In  (5.  12),  the  processes  §(uu,^,a)  and  r|(uu,H,a)  are  uncorrelated,  zero-mean  processes 
with  orthogonal  increments.  These  increments  have  variances  given  by 

2  2 

E[d5(uu,K,a)]  =  Eldr^uu.K.a)]  =  G(oj, n,a)  duu  dx  da  (5.13) 
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Heuristically,  we  can  continue  to  use  the  Fourier  transform  (5.  5)  to  represent  noise, 
if  we  interpret  X(uu,  x,ct)  and  Y(uu,  x,ct)  as  uncorrelated  random  functions  of  zero  mean 
satisfying 


E  X(u),K,a)  X^.K’.a')  =  E  Y(uu,x,a)  Y^’.x’.a') 

=  G(uu,K,a)  6(uu’ -  uu)  6(k' -  k)  6(a’ ~  a)  (5.14) 

If  we  wish  to  model  the  noise  as  a  superposition  of  surface -wave  modes  (uncorrelated), 
we  may  replace  G(uo,x,a)  by  a  sum  of  terms,  each  vanishing  off  the  corresponding  dis¬ 
persion  curve  of  the  mode  in  u u-x  space.  Since  these  separate  terms  in  G(uu,x,a)  will  not 
overlap,  the  representation  (5.  12)  (or  (5.  5)  )  may  be  written  as  a  sum  of  terms,  satis¬ 
fying  equations  like  (5.  13)  or  (5.  14)  for  each  mode. 

According  to  (5.  11),  any  assumed  power  spectral  density  in  uo-x  space  determines 
the  correlation  properties  of  the  noise  in  space  and  time.  For  example,  if  the  noise  power 
is  isotropic,  then  G(uu,  x,a)  =  Gq(uu,  x)/2tt,  is  independent  of  a  and  we  can  perform  the 
a- integration.  The  result  is 

00 

T(p,0,t)  =  JJ  Gq(uu,  x)  J0(np)  cos  out  duu  dx  (5.15) 

0 

where  J  is  the  Bessel  function  of  zero  order.  Note  that  the  correlation  function  is  also 
o 

isotropic.  If,  further,  this  isotropic  noise  is  propagating  in  a  single  surface-wave  mode, 
we  can  express  the  dispersion  curve  for  the  mode  in  the  form  x  =  x(uu).  Then  G()  will 
vanish  off  this  curve: 

G0(uu,  x)  =  P(uo)  6[x  -  x(tu)  ] 
and  we  obtain  the  single  integral 

CO 

T(p,t)  =  J  P(tu)  JQ[  px(uu)]  cos  out  duu  (5.16) 

0 

Since  Jo(0)  =  1,  we  see  that  p(uu)  is  the  one  “Sided  spectral  density  (per  radian)  of  the 
noise  at  a  given  location,  and  that  the  cross -spectral  density  for  the  noise  waveforms 
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at  two  points  a  distance  p  apart  is  completely  fixed  byP<uu)  and  the  assumptions  of 
isotropy  and  single -mode  dispersion 

More  generally,  we  may  expand  the  spectral  density  as  a  Fourier  series  in 
azimuth 


G(u),H,a>  =  ^  G0(ui,  h) 

03 

1  V  +  — 

+  —  [Gn  (u),h)  cos  na  +  Gn  (uush)  sin  na]  (5  17) 

n=l 

When  substituted  in  (5c  II),  we  can  perform  the  a-integration  to  obtain  the  correlation 
function  as  a  Fourier  series  in  0: 

oo 

r<p,e,t)  =  JJ  G0(u),  K)  J0(Kp)  COS  UJt  du)  dH 
0 

oo 

+  y;  cos  n  6  JJ  G  +(UU,H)  J  (Kp)  COS  <UUt 
n=l  0 

00 

■+•  'V  sin  n  6  fj  Gn  (od,h)  Jri(^p)  cos  (uut 
n~l  0 

Again,  die  double  integrals  may  be  reduced  to  single  integrals  (or  a  sum  of  single 
integrals)  by  the  assumption  of  a  single  mode  (or  a  sum  of  modes). 

The  processing  of  the  output  signals  of  an  array  can  be  discussed  in  terms  of 
processing  the  elastic  fie  d  itself.  For  example,  a  general  linear  functional  of  the 
field  has  the  form 


nrr 

—  )  duu  dH 


nn 


)  duu  dH 


(5^  18; 


JJJ  f(p,t)  uz(p,t)  dp  dt  (5. 19) 

—  oo 

for  some  fixed  function  f  An  actual  array  does  not  sample  the  field  continuously  in 
space,  but  discretely,  and  thus  is  described  by  an  f-function  which  is  a  sum  of  delta - 
functions  in  the  space  variable  However,  it  is  easier  to  see  what  is  happening  if  we 
use  space  integrals  instead  of  sums  over  sensors.  Formula  (5.  19)  describes  a  single 
number  based  on  the  field  for  all  values  of  time  Actual  processors  provide  a  function 
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of  real  time,  based  on  the  field  for  the  past  only  (i.e.  ,  they  use  realizable  filters).  If 
the  output  of  such  a  linear  processor  is  called  y(t),  we  may  write 


00 


00 


y(t)  =  JJ  do  J  ds  h(p,s)  Uz(p,t-s) 


(5.20) 


0 


where  the  function,  h(p,s),  completely  describes  the  processor.  In  practice,  many 
such  linear  processors  may  operate,  in  parallel,  on  the  same  input  field,  and  their 
outputs  may  be  combined  in  a  non-linear  manner.  However,  we  shall  analyze  only  the 
linear  aspect  of  the  processing.  For  an  array  of  sensors  at  points  p  =  pn, 
n  =  1,  2,  . .  .  ,  N,  we  have 


N 

h(Pjt)  =  Yj  hn(t)  o(p  -  pn) 


(5.21) 


n=l 


where  hp(t)  is  the  impulse  response  of  the  filter  applied  to  the  signal  at  the  n'th  sensor 
(including  the  response  of  the  seismometer  itself). 

In  order  to  express  the  processing  operation  (5.20)  in  ID--H  space,  we  use  the 
original  representation  (5.  1),  which  is  related  to  (5.  5)  by  means  of 


(5.22) 


Property  (5.  2)  is  accounted  for  by  definmg  X  to  be  even  in  uu,  and  Y  to  be  odd  in  id. 
Formula  (5.  1)  can  describe  noise,  as  well  as  signal,  if  we  employ  Eq.  (5.  14).  Now  if 
we  substitute  (5.  i)  in  (5.20),  we  find 


00 


(5.23) 


—  00 


where 


(5.24) 


—  00 


0 


Since  h(p,s)  is  real,  H(—  id,—  k)  -  H(uu,k),  and  (5.23)  maybe  written 
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2n  °° 

y(t)  -  2  Re  J  da  JJ  K  duu  dK  H(uu,H,a)  'Hou.K.a)  e 
0  0 


(5.25) 


Thus  H(aj,K,a),  for  values  of  ou  and  k  in  the  first  quadrant  of  the  oj-h  plane,  completely 
determines  the  processing.  For  the  array  mentioned  above,  we  have 


N 


H(oj,k)  =  £ 

n=l 


(5  26) 


where 


QO 


(5.27) 


is  the  frequency  response  of  the  filter,  hn(t). 

Formula  (5.25)  may  be  taken  as  the  starting  point  in  the  design  of  a  linear 
processor.  The  desired  response  in  uu-h  space  determines  H(ua,K,a),  which  then 
determines  h(p,t)  by  means  of  the  inverse  of  (5.24) 


dau  dK 


-  00 


0  0 


The  function  h(p,t)  is  then  approximated  by  a  discrete  sum  of  the  form  (5.21),  and  the 
resulting  processor  evaluated  for  performance,  aliasing  introduced  by  discrete  space¬ 
sampling,  etc.  Many  equivalent  procedures  also  suggest  themselves. 

As  an  illustration,  suppose  we  wish  to  sense  only  surface  waves  of  a  single  mode, 
described  by  a  dispersion  curve  h  -  h('jj),  in  a  given  frequency  band.  Then  we  may  put 
(for  oj  ^  0  and  k  >  0) 


(5.29) 
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This  will  accept  such  waves  isotropically,  with  a  frequency  weighting  function  F(uj). 
When  we  put  (5.29)  into  (5.28),  we  find 


h(p,t)  =  (2tt)”2  J  F(uo)  J0[pH(ui)]  e_li)t  daj  ,  (5.30) 

— 00 

where  F(—  uu)  =  F(uu)  .  The  processing  is  isotropic,  with  a  different  filter  for  each 
radial  distance.  For  example,  an  array  might  be  arranged  in  concentric  rings  of  sen¬ 
sors,  all  elements  of  a  ring  being  added  and  passed  through  a  filter  with  frequency 
response 


F(uu)  Jo[  PK  (uu)  ]  ,  (5.31) 

where  p  is  the  radius  of  the  ring.  As  a  special  case,  to  pick  out  a  particular  phase 
velocity,  c,  we  would  use  (5.31)  with  h(uu)  =  uu/c.  * 

If  we  wish  to  think  of  the  processing  in  terms  of  frequency  and  phase  velocity, 
we  simply  eliminate  H  by  a  change  of  variables: 

H  =  uu/c  (5.  32) 


where  0£c<“,  and  introduce 

2 

H  (uu,c,a)  =  ^  H(uj,  -,  a)  (5.33) 

i  c°  c 

Then  h(p,  t)  is  given  by 

2TT  ao 

h(p,t)  =  2  Re  (2tt)~3  J  da  JJ  duo  dc  H^uu.c  ,a)  e'l(uj/c)  COS  'a  0)  luJt  (5.34) 

0  0 

For  example,  if  we  make  independent  of  a  and  proportional  to  6(c  —  c^),  we  return  to 
(5.31)  with  h(uu)  =  uu/c0. 


*  Processing  schemes  which  accept  waves  within  a  fixed  range  of  phase  velocities  have 
been  designed  and  implemented  (Ref.  3).  They  have  proved  quite  useful  in  reducing 
microseismic  noi.se. 
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These  models  of  one  component  of  earth  motion  can  easily  be  extended  to  three- 

— *  ^ 

component  models  of  surface  motion.  The  displacement  is  now  a  vector,  u(p,t),  hence 
(5.  1)  is  generalized  to 


(5.35) 


where  the  transform,  <| J,  is  also  a  vector.  We  think  of  a  plane-layered  model  earth, 
and  draw  the  U)-H  diagram  for  the  underlying  half -space,  keeping  in  mind  that  (5.35) 
describes  the  motion  at  the  upper  (free)  surface  of  the  first  layer.  Using  the  results 
of  the  previous  sections,  we  know  that  eigenfunctions  are  possible  for  all  points  of  the 
continuum  portion  of  the  ud-h  space,  hence  \^(jl),k)  need  not  vanish  anywhere  in  this  por¬ 
tion.  For  fixed  u)  and  H,  the  component  of  ^  in  the  surface  and  perpendicular  to  H 
obviously  represents  the  amplitude  of  an  SH-wave,  and  we  know  that  there  is  just  one 
such  eigenfunction  (of  arbitrary  amplitude^  for  each  point  (ju,h).  The  other  two  com¬ 
ponents  of  i|f(uu,K)  describe  P—SV -waves,  and  we  know  that  there  are  two  independent 
P-SV  eigenfunctions  for  each  point  (uu,h).  Thus  we  may  choose  these  two  components 
of  l|j  arbitrarily,  and  we  have  shown  that  any  vector  amplitude  function,  t|/(uo,H)  is 
possible  for  the  continuum  portion  of  the  uu-H  space  In  the  surface -wave  portion  of  the 
uo -K  plane,  tjf  must  vanish  except  along  the  curves  describing  surface -wave  modes.  For 
each  such  mode,  the  polarization  is  a  given  function  of  jd  and  a,  hence  only  a  scalar 
function  of  frequency  can  be  assigned  arbitrarily  as  the  amplitude  of  each  mode. 

In  a  similar  fashion,  the  model  can  be  extended  in  depth  to  cover  the  general 
case,  by  writing 


(5.36) 


—  00 


For  each  eigenfunction,  the  z -dependence  is  fixed,  hence  it  is  not  appropriate  to  add  a 
fourth  variable  of  integration  conjugate  to  z,  and  no  new  degrees  of  freedom  are  present 
over  the  case  z  =  0.  Thus  a  statistical  specification  of  the  mode  structure  of  the  noise 
automatically  determines  the  correlation  properties  in  depth,  and  polarization,  as  well 
as  in  surface  location  and  time,  according  to  an  assumed  model  of  the  earth. 
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